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Abstract 


Given the limited underground water resources in Iran, careful calculation, proper use, 
regulation, and maintenance of these resources are crucial. One effective approach for 
managing and exploiting these resources optimally, both now and in the future, is through 
modeling. This study focuses on the challenges faced by groundwater in the Izadakhat 
Basin, employing GMS software and MODFLOW code to model and forecast water levels 
under both steady and unsteady conditions. Three scenarios are considered: continuation 
of current trends, a 20% decrease in precipitation, and a 20% increase in precipitation. The 
results indicate that, under current conditions, the water table is projected to drop by 2.786 
meters. A 20% decrease in precipitation would increase this decline to 3.77 meters, while 
a 20% increase in precipitation would reduce the drop to 1.77 meters. In other words, the 
water level fluctuates by about 1 meter with a 20% change in precipitation. Over the three- 
year period studied (2015-2017), the first scenario suggests a steady decline of 2.786 
meters every three years if current conditions persist. These findings highlight the 
importance of surface water storage in the region to replenish groundwater levels. 
Moreover, the analysis underscores the significant impact of human activity on both the 
quantity and quality of groundwater resources in the study area, emphasizing the need for 
improved resource management. 
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حکیده 


با توجه به محدودیت منابع آب زیرزمینی در ایران. محاسبه دقیق, استفاده صحیح» تنظیم و نگهداری این منابع از 
اهمیت بالایی برخوردار است. یکی از روش‌های Se‏ برای مدیریت و بهره‌برداری بهینه از این منابع در حال و آینده 
استفاده از مدل‌سازی است. این مطالعه بر چالش‌های مرتبط با منابع آب زیرزمینی در حوضه ایزدخواست تمرکز دارد 
و با استفاده از نرم‌افزار GMS‏ و کد MODFLOW‏ به مدل‌سازی و پیش‌بینی سطح آب در شرایط پایدار و IBLE‏ 
می‌پردازد. سه سناریو موردبررسی قرار گرفته است: ادامه روند فعلی. کاهش ۲۰ درصدی بارش و افزايش ۲۰ درصدی 
بارش. نتایج OUS‏ می‌دهد در شرایط فعلی. سطح آب زیرزمینی به میزان Y VAT‏ متر کاهش خواهد یافت. با کاهش ۲۰ 
درصدی بارش» این افت به ۳۰۷۷ متر خواهد رسید و با افزایش ۲۳ درصدی بارش این افت به ۱۸۷۷ متر کاهش 
می‌یابد؛ به‌عبارت‌دیگر: سطح آب با تغییر ۲۰ درصدی بارش تقریباً ۱ متر نوسان می‌کند. طی دوره سه‌ساله مطالعه‌شده 
COTATA E)‏ سناریوی اول OUS‏ می‌دهد در صورت ادامه شرایط فعلی» سطح آب هر سه سال یکبار به‌طور ثابت 
۲ متر کاهش خواهد یافت. این یافته‌ها بر اهمیت ذخیره آب‌های سطحی در منطقه برای تقویت سطح آب‌های 
زیرزمینی تأکید دارند. همچنین» تحلیل نتایج OUS‏ می‌دهد یکی از مهم‌ترین عوامل تغییر در کیفیت و کمیت منابع آب 
زیرزمینی در منطقه مطالعه‌شده عوامل انسانی است که باید در مدیریت این منابع موردتوجه قرار گیرد. 

کلیدواژه‌ها: نوسان سطح آب زیرزمینی» مدیریت pU‏ آب زیرزمینی. آبخوان. MODFLOWi.s‏ دشت 
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مقدمه 


قسمت اعظم آب شیرین موجود در کره زمین در زیر زمین قرار دارد. افزایش جمعیت. افزایش تقاضا بر 
منابع آب زیرزمینی را سبب شده است. به علت بهره‌برداری روزافزون از منابع آب زیرزمینی» مسائل جدیدی 
مطرح گردیده و موضوع مدیریت بهره‌برداری اقتصادی و مجاز و صیانت منابع نیز موردتوجه قرار گرفته 
است. ازآنجاکه منابع آب زیرزمینی کشور ایران محدود و درخور توجه است. این کمبود ایجاب می‌کند که در 
محاسبه مقدار آب موجود و استفاده صحیح و تنظیم و نگهداشت منابع آب زیرزمینی دقت کامل به عمل آید 
تا بتوان به‌طور مستمر و حتی برای نسل‌های بعدی از این منابع استفاده کرد. در واقع مکانیسم‌ها و شرایط 
مختلف مؤثر بر منابع آب زیرزمینی» تغییراتی در حجم ذخیره و سطح آب سفره‌های آبی به وجود می‌آورد که 
جهت برنامه‌ریزی کنترل و مدیریت صحیح Job‏ تحت نظارت دقیق پیش‌بینی‌های به‌موقع قرار گیرند Kamali‏ 
-(Ardakani, 2019)‏ یکی از روش‌های مفید برای مدیریت و بهره‌برداری بهینه در زمان حال و آینده استفاده از 


به‌طورکلی مدل‌ها به سه گروه مدل‌های فیزیکی. مدل‌های ریاضی و مدل‌های آنالوگ تقسیم می‌شوند. مدل‌های 
ریاضی ممکن است به‌صورت معین پا احتمالاتی يا ترکیبی از هر دو باشند. این مدل به مدل‌های تحلیلی و 
cote‏ نیز طبقه‌بندی می‌شوند )2019 (Kamali Ardakani,‏ که متداول‌ترین مدل‌های مدل‌سازی سطح تراز آب 
زیرزمینی از این گروه می‌باشند. 

در همین راستا در سال‌های اخیر در نقاط مختلف Lo‏ جهت مدل‌سازی و پیش‌بینی سطح تراز آب‌های 


زیرزمینی از مدل‌های متفاوتی استفاده شده است & (Knotters & Bierkens, 2001: Mackay, Jackson,‏ 


Wang, 2014; Shakya et al., 2022; Tao et al., 2022; Chidepudi et al., 2023; Mirzaei & Nazemi, 2011; 
برجسته‌ترین‎ MODFLOW Jw اما‎ Afruzi, & Zare Abyaneh, 2017; Akhoni Pourhosseini, Asadi, 2017 ) 


مدل در روش‌های ارزیابی منابع آب با استفاده از مدل‌های عددی جریان آب زیرزمینی است & (McDonald‏ 
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MODFLOW‏ انجام شده است )2018 (Rahmawati et al., 2013; Shakoor et al.,‏ که جریان را به‌صورت 


سه‌بعدی و به روش تفاضل محدود مدل می‌نماید. 


مدل MODFLOW‏ توسط سازمان زمین‌شناسی ایالات‌متحده در سال ۱۹۸۶ توسعه (McDonald & Cal‏ 
Harbaugh, 1988)‏ و به‌سرعت در نقاط مختلف دنیا جهت مدل‌سازی سطح تراز آب زیرزمینی مورداستفاده 
قرار گرفت (Rayne et al., 2001; Wang et al., 2008; Kushwaha et al., 2009; Shakoor et al., 2018; Pande‏ 

et al., 2022; Katibeh & ۲1۵/621, 2005; Heydari & Jabbari, 2022; Masoumi et al., 2021)‏ در همین 
راستا برخی محققین در خصوص شبیه‌سازی جریان و انتقال آلودگی به آب‌های زیرزمینی از این مدل استفاده 
کردند )2015 (Jovanovic et al., 2009: Mondal & Singh, 2009; Parameswari, & Mudgal,‏ برخی نیز 
جهت مدل‌سازی فرونشست )2021 (Rashvand et al., 2019; Rajabi et al.,‏ بررسی رابطه متقابل آب‌های 


سطحی و زیرزمینی( 2010 (Yanxun et al, 2011; Mehl & Hill,‏ و پیش‌بینی سطح آب زیرزمینی در 


سناریوهای مختلف تغییر اقلیم )2022 (Patil et al, 2020: Abbasnouvinpour et al,‏ و خشکسالی 


(Atayizadeh & Chit Sazan, 2018; Porhaghi, et al. 2014)‏ از مدل مذ کور استفاده کردند و مورد تان 


قرار دادند. 


بنابراین در این راستا و با توجه به اهمیت این موضوع و مشکلاتی از قبیل اقلیم گرم و خشک. پتانسیل بالای 
تبخیر و تعرق و شوری بالای آب‌های زیرزمینی در استان فارس که تهدیدی جدی برای شور شدن خاک و 
ee‏ مق oae eal‏ کر که ی رف cul‏ فط لے امت adliuc Salts‏ در جوت اسان 
فارس» واقع در بخش مرکزی شهرستان زرین‌دشت انتخاب شده است. این منطقه یک حوضۀ بسته است که 
کمبود آب مطلوب و شور شدن آب‌وخاک در سال‌های اخیر مانع جدی برای کشاورزی و معیشت مردم شده 
است. به‌طوری‌که ویژگی‌های زمین‌شناسی و ژئومورفولوژی باعث تغییر کمیت و کیفیت منابع آب شده و 
مشکلات مضاعفی را برای این منطقه به وجود آورده است. همچنین در سال‌های اخیر با افزايش جمعیت» 
چاه‌هایی که برای مصارف خانگی استفاده می‌شوند. از حالت کم‌عمق به نیمه‌عمیق تبدیل و تعداد قابل‌توجهی 


از چاه‌های بهره‌برداری نیز خشک شده‌اند و به‌دلیل کمبود آب در چاه‌های کم‌عمی. کشت و زرع کاهش پیدا 
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کرده است. حتی در بعضی مناطق استفاده از آب شور جهت کشاورزی همراه با کیفیت پایین خاک و سایر 
شیوه‌های غلط مدیریت منجر به تأثیرات زیان‌آور به منابع خاک شده است. این در حالی است که با توجه به 
بحران. کمبود و شوری col‏ در منطقه اقدامات مناسبی جهت مدیریت منابع آب انجام نشده است؛ بنابراین با 
توجه به مشکلات do‏ کور. در پژوهش حاضر سعی بر این است که ضمن بررسی مشکلات ایجاد شده سطح 
تراز آب زیرزمینی حوضه ایزدخواست با استفاده از نرم‌افزار GMS‏ با کد MODFLOW‏ در شرایط ماندگار و 
غير ماندگار مدل‌سازی و در سه سناریوی تداوم شرایط جاری. کاهش و افزایش ۲۰ درصدی بارش پیش‌بینی 


)- مواد و روش‌ها 


-V-Y‏ منطقه موردمطالعه 


حوضه آبریز ایزدخواست با مختصات ۲۹۳ £47 ۵۳۳ تا YA"‏ ۸ ۵4۳ طول شرقی و YA? V OY"‏ تا ۲۰۳ YA‏ 
۸ عرض شمالی یکی از کفه‌های بستة حوضه آبریز رودخانۀ مند واقع در استان فارس است (شکل O‏ 
مساحت این حوضه ۱۳۷۱/۳ کیلومترمربع است که ارتفاعات ۸۷۹/۱ و دشت 1٩۱/۷‏ کیلومترمربع از کل 
مساحت حوضه را به خود اختصاص داده‌اند. حداکثر ارتفاع در حوضه ۲۱۸۲ و حداقل pur‏ ۱۰۲۹ متر 
است. این حوضه دارای ميانگین بارش سالانه ۲۲۲/۱ میلی‌متر میانگین سالانه دما ۱۹/۵ درجه و متوسط تبخیر 


سالانه از تشتک ۲۹۵۵/۹ میلی‌متر است(2 2021 (Ansari et al.,‏ . 


از نظر هیدرولوژی در این حوضه رودخانه دائمی وجود نداشته و AS‏ جریان‌ها موقت و فصلی می‌باشند که 
تمامی جریان‌های سطحی منطقه به‌صورت همگرا از حاشیه به سمت مرکز دشت = BL OL‏ و درنهایت در 
مرکز و غرب دشت تجمع می‌پابند. در مرکز دشت آب‌های راکد پس از تبخیر» اراضی شور و غیرقابل استفاده 
را نمایان می‌سازند. عمده‌ترین مسیل‌های منطقه مسیل خشک آنگبینه واقع در غرب روستای دوبران» مسیل 
خشک چاه زبر و مسیل خشک گلکویه (گلکان) هستند (Ansari et al, 2021 b)‏ . این جریان‌ها بر روی 


سازندهای مقاوم‌تر مانند آسماری جهرم که دارای شکستگی‌ها ودرز و شکاف فراوان هستند بستر خود را به 
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شکل دره‌های ژرف ۷ مانند تبدیل ساخته‌اند اما بر روی سازندهای سست‌تری مانند آغاجاری دارای بستری 


پهن هستند و اغلب آن‌ها پس از ورود به دشت بستر مشخصی ندارند. 


-7Y‏ روش انجام پژوهش 
جهت مدل‌سازی سطح تراز آب‌های زیرزمینی در این تحقیق از نرم‌افزار GMS kl;‏ و MODFLOW us‏ 
استفاده می‌شود. این نرم‌افزار که مدل l5 MODFLOW‏ 22 2 حای داده است» با امکانات مختلفی که به ol‏ 


افزوده شده است نسبت به خود نرم‌افزار MODFLOW‏ توانایی بیشتری در وارد کردن des‏ ساخت و 


تحلیل مدل دارد. 


چاه های مشاهداتی e‏ 
نقاط ارتفاعی حوضه ایزدخواست 
Val‏ 


uc 
- High : 2182 


ات سس سح 1029 : Low‏ 
26 195 13 3/2565 0 


شکل ۱- موقعیت منطقه موردمطالعه 


Fig.1. Location of the study area 
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کد MODFLOW‏ بر ol‏ حل معادلات > ES‏ آب‌های زیرزمینی Cul‏ به‌طوری‌که حرکت سه‌بعدی آب 
زیرزمینی با چگالی ثابت به‌وسیله معادله دیفرانسیلی جزئی (رابطه ۱) با استفاده از روش تفاضل محدود و بر 


اساس معادله پیوستگی حل می‌شود. 
on] a ðh a an] ô.,‏ 
رابطة [vy = + az [kzz = -W Ss — = d [xx = ay"‏ 


تغذیه و تخلیه o >> 5 wl‏ مخصوص مواد متخلخل. # زمان (روز) 9 X‏ ۷و 2 بیانگر مختصات است. 


برای ایجاد مدل جریان در نرم‌افزار GMS‏ از روش مدل مفهومی | استفاده شده است. در این روش ابزارهای 
GIS‏ و مدل Map‏ به‌کار برده می‌شوند. تمامی داده‌های موجود در منطقه ازجمله موقعیت منابع تغذیه و تخلیه 
et db‏ اها و d‏ خی بارا aula‏ لباز قا cultas‏ هدور فر قان سل و SU‏ 
دیگر موردنیاز برای مدل‌سازی در سطح مدل مفهومی موردبررسی قرار گرفت و به نرم‌افزار داده شد. سپس 
شبکه ایجاد و مدل مفهومی به مدل شبکه‌ای " تبدیل می‌گردد. محدوده دشت به مدل وارد می‌شود. سپس هر 
یک از ویژگی‌های آبخوان ale‏ تغذیه و تخلیه و چاه‌های مشاهده‌ای در یک پوشش جداگانه تعریف 
می‌شوند. در مدل مفهومی محدوده موردمطالعه» موقعیت چاه‌های بهره‌برداری و مشاهده‌ای به مدل داده شد. 
به‌طورکلی ساختار مدل مفهومی آبخوان شامل محدوده مدل‌سازی و توزیع adsl‏ پارامترهای هیدروژئولوژیک؛ 
adas‏ چاه‌های بهره‌برداری و میزان آب برگشتی آن‌ها. چاه‌های مشاهداتی. میزان تغذیه از سطح به آبخوان و 
شرایط مرزی آبخوان خواهد بود. 

در ادامه» اطلاعات جمع آوری شده در مدل مفهومی شامل مقادیر پارامترهای فیزیکی (ضخامت لایه آبدان 


ضریب تخلخل و غیره) و هیدرودینامیکی (ضریب هدایت هیدرولیکی و قابلیت انتقال) در آبخوان و هم 


چنین مقادیر تغذیه (نفوذ ناشی از بارش و Ol‏ برگشتی ناشی از مصارف) و تخلیه از آبخوان (پمیاژ چاه‌ها) با 


' Conceptual Model Approach 
2. Grid Model 
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استفاده از داده‌های سازمان سهامی آب منطقه‌ای استان فارس. سازمان هواشناسی و DEM‏ ۵ متر منطقه به 
مدل ریاضی منتقل شد و هرکدام توسط نماد مختص به خود نمایش داده شد. 

با توجه به نبود اطلاعات هدایت هیدرولیکی در منطقه موردمطالعه داده‌های مربوط به ضریب قابلیت انتقال 
محدوده موردنظر از سازمان سهامی آب منطقه‌ای استان فارس تهیه گردید و از تقسیم ضریب قابلیت انتقال به 
ضخامت a‏ اشباع آبخوان. هدایت هیدرولیکی قسمت‌های مختلف دشت تخمین زده شد. مراحل آماده‌سازی 
پارامترهای هیدرودینامیکی در نرم‌افزار ArcGIS‏ انجام شد. همچنین L‏ توجه به مطالعات انجام‌شده توسط 
وزارت y‏ ضصریب ذخیره آبخوان Cans‏ ایزدخحواست به‌طور متوسط 0 درصد در نظر گرفته شد و به‌مدل 
معرفی گردید. از طرفی Lb‏ توجه داشت که میزان ضریب ذخیره و آبدهی ویژه یکی هستند و به یک معنا 


در مرحلۀ بعد با توجه به مقدار بارش و ضریب نفوذپذیری مقدار تغذیه ناشی از بارندگی مشخص شد. پس 
از ol‏ چاه‌های کشاورزی دشت مشخص گردید. با توجه به ایستگاه‌های باران‌سنجی داخل و پیرامون دشت 
برگرفته از سایت سازمان هواشناسی استان فارس, مقدار متوسط بارندگی استخراج شد و با توجه به مطالعات 
ah‏ انجام شده ضریب نفوذ بارندگی در تمام سطح دشت ply‏ ۵ درصد در نظر گرفته شد (سازمان سهامی 
آب منطقه‌ای استان فارس,۱۳۹۶) (شکل ۲). لازم به ذکر است که مدل‌سازی سطح تراز آب زیرزمینی فقط در 


دشت منطقه موردمطالعه انجام شده است. 


در dle ye‏ بعد چاه‌های مشاهداتی به‌وسیله نرم‌افزارهای GIS‏ و Excel‏ آماده شدند و بعد از تبدیل به فرمت 
(Shape Files) L (Text Files)‏ به GMS‏ در سطح مدل مفهومی (در ماژول نقشه) وارد و پس از آن» به 
مدل شبکه‌ای سه‌بعدی MODFLOW‏ تبدیل شدند. سطح آب به نقاط مشاهده‌ای (شکل ۳) که از قبل تعیین 
شده بودند. تخصیص داده شد. سپس تعداد ۱۲۳ حلقه چاه بهره‌برداری (شکل (E‏ در محدوده دشت 
ایزدخواست با استفاده از اطلاعات سازمان سهامی آب منطقه‌ای استان فارس با حجم کل تخلیه سالانه از 


جاه‌های موجود در این آماربرداری حدود ۰ میلیون مترمکعب محاسبه گردید. 
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شکل ۲- تغذیه از سطح آبخوان 


Fig.2. Recharge from the aquifer surface 


شکل ۳- پراکنش چاه‌های مشاهده‌ای در منطقه موردمطالعه 


Fig.3. Distribution of observation wells in the study area 
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شکل ۶- پراکنش چاه‌های بهره‌برداری در منطقه موردمطالعه 


Fig.4. Distribution of extraction wells in the study area 
استخراج‎ GIS و با استفاده از نرم‌افزار‎ DEM در مرحلة بعد اطلاعات رقومی تویوگرافی سطح زمین از داده‌های‎ 
به مدل وارد شد. سپس با استفاده از درون‌یابی» کلیه نقاط به‌طور مستقیم به‎ shp گردید و به‌صورت فایل‎ 
شبکه مدل نسبت داده شد. نقشه مدل رقومی ارتفاعی و سنگ کف آبخوان نیز به کمک داده‌های ژئوفیزیک و‎ 


نقشه تغییرات ضخامت آبخوان ساخته شد (شکل‌های ۵ و .)١‏ 


پس از مشخص کردن پارامترهای زمانی و پارامترهای هیدرولوژیک و هیدروژئولوژیک برای هر سلول مدل. 
توصیف دقیق شرایط مرزی با استفاده از مناسب‌ترین بسته‌های MODFLOW‏ صورت می گیرد. در این مرحله از 
تحقیق» شرایط مرزی سفره عمدتاً طبق نتایج مطالعات پایه آب‌های زیرزمینی و بر اساس نقشه‌های تراز آب» 
داده‌های مشاهداتی دشت و شرایط زمین‌شناسی و هیدروژئولوژیکی سفره تعیین شده است. مناطقی که 
هم‌جوار سازندهای نفوذپذیر قرار داشتند به‌عنوان مرزهای نفوذناپذیر تعیین شدند. سال ul‏ | ۱۳۹۵-۱۳۹۶ به 
علت کامل بودن اطلاعات در این سال به‌عنوان سال ab‏ شبیه‌سازی در نظر گرفته شدند و مقادیر ارتفاع سطح 
ایستابی در مهرماه AE‏ که ماه حداقل سطح آب در منطقه بود به‌عنوان شرایط اولیه در حالت پایدار انتخاب 
گردیدند. همچنین شبیه‌سازی سطح آب و واسنجی مدل در حالت غیر ماندگار در مهرماه سال آبی -WA^‏ 


۱۶۰۳ تابستان‎ Y شماره‎ OY نشریه جغرافیا و مخاطرات محبطی» جلد‎ \VY 


بنابراین پس از تعیین تمام داده‌های موردنیاز کد MODFLOW-2000‏ مدل آماده اجرا است. در ابتدا با استفاده 
از Model Checker‏ خطاهای احتمالی بررسی می‌شود و تصحیح می گردد. در مرحله بعد یکی از روش‌های 
محاسباتی MODFLOW‏ را انتخاب و مدل اجرا می‌شود که در این مطالعه از روش تکرار ۳062 استفاده 
گردید. معمولاً در اولین اجراء برازش مناسبی بین بار هیدرولیکی محاسبه‌شده و اندازه‌گیری‌شده به دست 
نخواهد آمد. بنابراین برای رسیدن به یک برازش قابل‌قبول Lb‏ خطاهای احتمالی موجود در فرایند مدل‌سازی 
شناسایی و تا حد ممکن رفع گردد و مدل واسنجی شود که در این مطالعه از روش واسنجی دستی و همچنین 


از کد PEST‏ جهت واسنجی اتوماتیک مدل استفاده گردید. 


شکل ۵- نقشه DEM‏ منطقه موردمطالعه 


Fig.5. DEM map of the study area 
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Fig.6. Map of the aquifer bottom in the study area 
واسنجی مدل را می‌توان برای شرایط ماندگار و غیرماندگار و یا هر دو انجام داد. برای ارزیابی نتایج واسنجی‎ 
جهت محاسبه اختلاف سطح ایستابی شبیه‌سازی‌شده و سطح‎ RMSE از پارامتر‎ GMS از امکانات نرم‌افزار‎ 
۱۳۹۶ سال‎ ole ge ایستابی مشاهده‌شده در چاه‌های پیزومتری و تجزیه‌وتحلیل آن‌ها استفاده شد. در این تحقیق‎ 
ماه پایه برای واسنجی مدل در حالت رژیم پایدار انتخاب شد. در این حالت بعد از اجرای مدل» ابتدا‎ ol pew 
هدایت هیدرولیکی ميزان تغذیه و تخلیه استفاده‎ polis جهت واسنجی‎ (PEST) از روش واسنجی اتوماتیک‎ 
شد و محدوده تقریبی پارامترها مشخص شد و برای رسیدن به واسنجی دقیق‌تر از روش سعی و خطا (دستی)‎ 
شده ارتفاع سنگ بستر و مرزهای دشت از نظر هد هیدرولیکی ثابت‎ Sd استفاده شد که علاوه بر پارامترهای‎ 
سطح ایستایی اندازه‌گیری شده توسط مدل با مقادیر اندازه‌گیری شده در‎ calis نیز تغییر داده شد تا‎ 


چاه‌های مشاهده‌ای بسیار نزدیک گردد. 


مشخص شد. در این مرحله به واسنجی ضریب ذخیره و ميزان تغذیه و تخلیه پرداخته شد که عوامل تغذیه و 


تخلیه در یک دوره Y‏ ساله ۱۳۹۶-۱۳۹۲ با ۱۲ دوره تنش به مدل وارد گردید و به روش آزمون سعی و خطا 


مورد واسنجی قرار گرفت و همچنین از روش واسنجی اتوماتیک 1 برای تخمین Ol‏ ضریب ذخیره 
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استفاده گردید. در خلال این فرآیند shee‏ واسنجی سطح ایستابی مشاهده‌شده در چاه‌های مشاهده‌ای برای ۱۲ 
eb‏ دوره شبیه‌سازی بود. همچنین نتایج حاصل از این مرحله با کمک پارامتر RMSE‏ موردبررسی قرار 
گرفت. درنهایت در مدل واسنجی شده آبخوان محدوده موردمطالعه. با تعیین عدم قطعیت داده‌های $2329« 
پارامترهای مدل‌سازی در دامنه‌ی معقول و مشخص تغییر داده شدند. حساس‌ترین پارامتر پارامتری است که 
کمترین تغییر در آن موجب بیش‌ترین تغییر در واسنجی شود. در این مطالعه با استفاده از PEST az.‏ از 


نرم‌افزار GMS‏ حساسیت مدل نسبت به پارامترهای مختلف سنجیده شد. 


جهت اجرای سناریوهای پیش‌بینی چون خشکسالی و کمبود بارندگی بر منابع آب زیرزمینی و رطوبت خاک 
و جریان رودخانه موثر است )2014 Porhaghi, et al.‏ و با توجه به اين‌که با پیش‌بینی و مدل‌سازی منابع آب 
زیرزمینی می‌توان نقاط ple‏ به کاهش و افزايش تغذیه را شناسایی و برای مدیریت بهتر از آن‌ها استفاده 
نموده بنابراین از آمار ۱۳ JUS‏ بارندگی موجود در منطقه که از ادار کل هواشناسی استان فارس گرفته hd‏ 
استفاده گردید. بر Gb‏ این آمار میانگین بارندگی حوضه از سال ۱۳۸۷ تا سال ۰۱۳۹۹ ۲۱۷/۵۷ میلی‌متر است. 
کم‌ترین میزان بارندگی در حوضه مربوط به سال ۱۳۸۷ با ۵۰/۱۲ میلی‌متر بارش و بیش‌ترین میزان بارش 
مربوط به سال ۱۳۹۵ با ٤۹۹/٤‏ میلی‌متر بارش گزارش شده است. پس با توجه به میانگین و حداقل بارش در 
این دور ۱۳ ساله سناریوهای پیش‌بینی با تداوم شرایط کنونی و کاهش و افزايش ۲۰ درصدی بارش به مدل 


اعمال شد. 
qs -۳‏ و بحث 


با توجه به نتایج مقادیر واسنجی شدۀ هدایت هیدرولیکی به دست آمده پس از کالیبراسیون مدل در محدوده 
مدل‌سازی و در lei‏ پایدار OLS‏ می‌دهد که بیش‌ترین مقادیر در قسمت مرکزی حوضه متمرکز شده (شکل 
۷ و بیش‌ترین میزان AS‏ نیز در مرکز و جنوب دشت قرار دارند (شکل ۸) که با توجه به شکل و 
توپوگرافی حوضه مقادیر قابل قبولی به‌حساب می‌آیند؛ بنابراین در مرحلۀ بعد مقایسه بین سطح آب مشاهده 
شده و محاسبه شده در محل چاه‌های مشاهداتی در OLL‏ دوره واسنجی پایدار انجام گرفت و تفاوت آن‌ها 


مشخص شد (جدول ۱). همچنین احتلاف سطح آب در چاه‌های مشاهده‌ای با محاسباتی در دوره واسنجی 
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پایدار به‌صورت نمودار میله‌ای ترسیم شده است (شکل .)٩‏ با توجه به جدول ۱ و شکل ٩‏ می‌توان دریافت 
که اختلاف تراز چاه‌های مشاهداتی و محاسباتی در تمام چاه‌ها کم‌تر از یک متر است به‌استثنای چاه شمارۀ ۷ 
که این اختلاف ۱/۰۵ متر است که نتایج قابل قبولی می‌باشند که می‌توان این مدل‌سازی را در شکل ۱۰ 
مشاهده نمود. همان‌طور که گفته شد چون حداکثر اختلاف بین palis‏ مشاهداتی و محاسباتی تراز آب» 1۱ 
در نظر گرفته شده است در این شکل رنگ سبز پیزومترها OUS‏ می‌دهد که معیار واسنجی در حالت ماندگار 


cule, ELT‏ شده و تفاوت مقادیر مشاهداتی و محاسباتی کمتر از (EV)‏ است. 


شکل ۷- مقادیر واسنجی‌شده هدایت هیدرولیکی در محدوده مدل‌سازی طی کالیبراسیون مدل در شرایط پایدار 


Fig.7. Calibrated hydraulic conductivity values in the modeling area during steady-state calibration 


شکل -A‏ مقادیر واسنجی شد: تغذیه طی کالیبراسیون مدل در شرایط پایدار 


Fig.8. Calibrated recharge values during steady-state model calibration 
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حدول -Y‏ تراز مشاهداتی و محاسباتی در حالت پایدار 


Tablel- Observational and computational balance in steady state 


چاه مشاهداتی (متر) محاسباتی (متر) 
4VA4/AA ava W1‏ 
44Y/of 44/۱۲ W2‏ 
CASTA w3‏ 4/۲ 
AVY/AV 4/1۸ W4‏ 
REV W5‏ ۱۱۳/۲۱ 
W6‏ ۳۷۱ 1/11 
W7‏ 134/1۷ 11۲ 
W8‏ ۱ 100/۳ 
AL ۹۸/۸ W9‏ 
avy W10‏ 47771 
1 3 
0.5 7 
olm | | ۱ ۱ " | 8‏ 4 
9 
E -1.5‏ 
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دوره واسنجی پایدار‎ ULL شکل 4- مقایسه بین سطح آب مشاهداتی و محاسباتی در‎ 


Fig.9. Comparison between observed and calculated water levels at the end of the steady-state 
calibration period 
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شکل ۱۰- جاه‌های مشاهداتی آبخوان در محدوده مدل‌سازی 


Fig.10. Observation wells of the aquifer in the modeling area 
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شکل ۱۱- مقادیر مدلسازی و مشاهداتی آبخوان در محدوده مدل‌سازی 
Fig.11. Modeled and observed values of the aquifer in the modeling area‏ 
در شرایط غیرماندگار نیز مانند شرایط ماندگار واسنجی به‌وسیله ترکیبی از دو روش خودکار PEST‏ و سعی و 
خطای دستی. انجام شده و مقادیر آبدهی ویژه و هدایت هیدرولیکی برای سطح محدوده مدل‌سازی واسنجی 
شد. همچنین در حالت غیرماندگار حداکثر اختلاف بین polis‏ مشاهداتی و محاسباتی تراز آب EN)‏ در نظر 


گرفته شد (شکل QY‏ 
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بر اساس شکل ۳ می‌توان مشاهده نمود که در حالت غیرماند JS‏ نیز در تمام جاه‌ها معیار واسنجی ELT‏ 
cule,‏ شده و تفاوت مقادیر مشاهداتی و محاسباتی کمتر از (EV)‏ است. پس با دقت زیادی می‌توان به نتایج 


OY‏ و سناریوی سوم با افزايش ۲۰ درصدی بارش (شکل (VE‏ در نظر گرفته شد. 


scenario-1 


= 
روند کنونی‎ Ub شکل ۱۲- سناریوی پیش‌بینی سطح تراز آب زیرزمینی با‎ 
Fig.12. Scenario predicting groundwater level with the continuation of the current trend 


scenario-2 
22.880 


i 
شکل ۱۳- سناریوی پیش‌بینی سطح تراز آب زیرزمینی با کاهش ۲۰ درصدی بارش‎ 


Fig.13. Scenario predicting groundwater level with a 20% reduction in precipitation 
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scenario-3 


909.293 


L 
سناریوی پیش‌بینی سطح تراز آب زیرزمینی با افزایش ۲۰ درصدی بارش‎ M شکل‎ 
Fig.14. Scenario predicting groundwater level with a 20% increase in precipitation 
تغییر میزان بارش در منطقة‎ UL درنهایت نتایج پیش‌بینی شرایط در سه مورد سناریو به دست آمد که بر‎ 
موردمطالعه تعریف شده است. در سناریوی اول امتداد شرایط حاضر و در سناریوهای دوم و سوم میزان‎ 
کاهش و افزایش ۲۰ درصدی بارش در نظر گرفته شد؛ بنابراین با توجه به داده‌های تراز چاه‌های مشاهداتی‎ 
و حداکثر ۱۱۲۳/۶۵ را نشان می‌دهند در نظر گرفتن امتداد شرایط حاضر سطح‎ ٩۱۱/۰۲ سال ۱۳۹۹ که حداقل‎ 
۲۰ می‌دهد. با حالت کاهش‎ OUS می‌رسد که ۲/۷۸۲ متر کاهش را‎ ٩۰۸/۲۷۶ آب در حالت حداقل به مقدار‎ 
متر و با افزایش ۲۰ درصدی بارش این اختلاف به ۱/۷۷ متر می‌رسد؛‎ Y ۸۷۷ درصدی بارش این اختلاف به‎ 
یعنی با کاهش و افزایش بارش تقریباً ۱ متر بر میزان سطح تراز آب در سناریوی امتداد شرایط کاسته و‎ 
UL. افزوده خواهد. میزان حداکثر سطح تراز نیز این شرایط را دارا می‌باشد که با توجه په در نظر گرفتن سه‎ 
۲/۷۸۲ سطح تراز (۱۳۹۹-۱۳۹۶) سناریوی اول گویای این مطلب است که با امتداد شرایط کنونی هر سه سال‎ 
چاه‌های‎ DL ۱۰ متر از میزان سطح تراز آب کاسته خواهد شد. صحت این موضوع را می‌توان در گزارش‎ 
منطقه‌ای استان فارس (۱۳۹۶) مشاهده کرد. بر اساس گزارش‎ Ol بهره‌برداری مربوط به سازمان سهامی‎ 
چاه بهره‌برداری موجود در منطقه» ۶۳ چاه متروکه و خشک می‌باشند که‎ VIT مذکور در سال ۱۳۸۸ از تعداد‎ 
مشاهداتی به ۱۶۵ چاه متروکه و خحشک رسیده است که از پایین‎ ole ۳۰۲ این رقم در سال ۱۳۹۸ از تعداد‎ 


رفتن سطح تراز در ده سال مذکور حکایت دارد. پر این اساس می‌توان همجون پژوهش‌های قبل قدرت مدل 
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۷ در مدل‌سازی و پیش‌بینی سطح تراز آب زیرزمینی در سناریوهای مختلف را تأیید نمود 
(Rahmawati et al., 2013; Shakoor et al., 2018; Porhaghi et al., 2014; Heydari & Jabbari, 2022 )‏ 

نکتۀ قابلتوجه دیگر نیز در سناریوها وجود حداقل سطح تراز آب‌های زیرزمینی در شرق و غرب حوضه 
می‌باشد که می‌تواند آسیب‌پذیری منابع آب در این قسمت‌ها را بیشتر نماید؛ بنابراین براساس این نتایج ذخيرة 
آب‌های سطحی در منطقه جهت افزايش سطح تراز آب‌های زیرزمینی ضروری است. همانند پژوهش 
45(Shakoor et al., 2018)‏ تغذیه مصنوعی جهت افزایش سطح تراز را پيشنهاد دادند؛ بنابراین اولین اقدام در 
جهت iue‏ منابع آب سطحی آگاهی از وضعیت کیفی منابع آب سطحی و شناسایی منشأً شوری منابع 


آب‌وخاک می‌باشد. 


بر این اساس از مهم‌ترین عوامل تغیر کیفیت و کمیت منابع آب زیرزمینی در منطقه موردمطالعه عوامل انسانی 
می‌باشد. به‌طوری که می‌توان به‌طور واضحی SU‏ فعالیت انسانی را در گزارش‌های سازمان سهامی آب منطقه- 
ای (۱۳۹۶) در این حوضه مشاهده کرد. بر طبق گزارش این سازمان در دور مطالعاتی (۱۳۹۸-۱۳۸۸) هرساله 
نزدیک به ۳۰ میلیون مترمکعب آب از چاه‌ها استخراج می‌شده است که به‌مرور به‌دلیل حشک شدن برخی 
چاه‌ها این مقدار به ۲۰ میلیون مترمکعب در سال ۱۳۹۸ رسیده است که با توجه به بیلان منفی منابع آب 
زیرزمینی در این حوضه سالانه اضافه برداشت از مخازن col‏ به ميزان ۰/۱۷۱ میلیون مترمکعب وجود دارد که 
این موضوع به‌تدریج از سطح تراز آب زیرزمینی در این مخازن کاسته و بر میزان هدایت الکتریکی آب افزوده 
شده است. )2022 (Heydari & Jabbari,‏ نیز در پژوهشی که در دشت مرودشت استان فارس انجام دادند. به 
این نتیجه رسیدند که میزان برداشت بی‌رویه آب زیرزمینی Mad)‏ مستقیمی با افت آب‌های زیرزمینی دارد؛ که 
این امر می‌تواند تو b uS‏ شوری منابع col‏ خاک منطقه کمک (Ansari et al, 2023; Ansari £S‏ 
s&jabbari,2023)‏ و ازآنجاکه این آب‌ها در محدودة زمین‌های زراعی استفاده شده است و در جنوب شهر 
دبیران بر افزایش شوری خاک زمین‌های زراعی افزوده شده است می‌توان با دقت بیشتری به این موضوع 


توجه نمود. ولی براساس کاهش سطح تراز و افزایش هدایت الکتریکی از سال ۱۳۹۲ تا ۱۳۹۹ می‌توان تأثیر 
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فعالیت‌های انسانی را در تغییر کیفیت منابع آب زیرزمینی YEN.‏ درصد و تأثیر عوامل طبیعی را 70/۳ درصد 


(Ansari et al., 2023) تخمین زد‎ 


*- نتیجه گیری 


نتایج به‌دست آمده در این پژوهش OLE‏ می‌دهد که با در نظر گرفتن تداوم شرایط حاضر سطح آب ۲/۷۸ 
مقر کاهش می‌یابد. با حالت کاهش ۲۰ درصدی بارش این احتلاف به ۸۷۷ Y‏ متر و با افزایش ۲۰ درصدی 
بارش این اختلاف به ۱/۷۷ متر می‌رسد؛ یعنی با کاهش و افزایش بارش تقریباً ۱ متر بر میزان سطح تراز آب 
در سناریوی تداوم شرایط کاسته و افزوده خواهد. میزان حداکثر سطح تراز نیز این شرایط را دارا می‌باشد که 
با توجه به در نظر گرفتن سه AUS‏ سطح تراز (۱۳۹۲-۱۳۹۶) سناریوی اول گویای این مطلب است که با 
تداوم شرایط کنونی هر سه سال ۲/۷۸۰ متر از میزان سطح تراز آب کاسته خواهد شد و اگر بارش ۲۰ درصد 
کاهش LL‏ این رقم به ۳/۷۷ خواهد رسید و در شرایط بسیار خوش‌بینانه (سناریوی سوم) هر سه سال ۱/۷۷ 
متر از سطح تراز آب زیرزمینی منطقه کاسته خواهد شد که بسیار برای حوضه چالش‌برانگیز می‌باشد و 


مدیریت مناسبی را می طلبد. 


بنابراین براساس این نتایج ذخیرة آب‌های سطحی در منطقه جهت افزایش سطح تراز آب‌های زیرزمینی 
ضروری است. همچنین تحلیل نتایج نشان داد که از مهم‌ترین عوامل تغیر کیفیت و کمیت منابع آب زیرزمینی 
در منطقه موردمطالعه عوامل انسانی مانند افزایش زمین‌های زیرکشت و به دنبال آن افزایش چاه‌های بهره- 
برداری و افزایش برداشت آب زیرزمینی می‌باشد که جهت مدیریت این منابع باید مدنظر قرار گیرد. 

درنهایت نتایج نشان idas‏ قدرت مدل MODFLOW‏ در مدل‌سازی و پیش‌بینی سطح تراز آب زیرزمینی در 
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